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We study the equation of state for symmetric nuclear matter using a ring-diagram approacti in 
which the particle-particle hole-hole (pphh) ring diagrams within a momentum model space of dec- 
imation scale A are summed to all orders. The calculation is carried out using the renormalized 
low-momentum nucleon-nucleon (NN) interaction Viow-k, which is obtained from a bare NN po- 
tential by integrating out the high-momentum components beyond A. The bare NN potentials of 
CD-Bonn, Nijmegen and Idaho have been employed. The choice of A and its influence on the single 
particle spectrum are discussed. Ring-diagram correlations at intermediate momenta (A: ~ 2 fm~^) 
are found to be particularly important for nuclear saturation, suggesting the necessity of using a 
sufficiently large decimation scale so that the above momentum region is not integrated out. Using 
Viow-k with A ~ 3 fm~^, we perform a ring-diagram computation with the above potentials, which 
all yield saturation energies E/A and Fermi momenta A:^' considerably larger than the empirical 
values. On the other hand, similar computations with the medium-dependent Brown-Rho scaled 
NN potentials give satisfactory results of E/A ~ — 15 MeV and fc^' ~ 1.4 fm~^. The effect of this 
medium dependence is well reproduced by an empirical 3-body force of the Skyrme type. 

PACS numbers: pacs 



I. INTRODUCTION 



Obtaining the energy per nucleon (E/A) as a func- 
tion of the Fermi momentum (kp) for symmetric nuclear 
matter is one of the most important problems in nu- 
clear physics. Empirically, nuclear matter saturates at 
E/A ~ —16 MeV and kp — 1-36 fm~^. A great amount 
of effort has been put into computing the above quanti- 
ties starting from a microscopic many-body theory. For 
many years, the Brueckner-Hartree-Fock (BHF) theory 
[Hi Si S| was the primary framework for nuclear matter 
calculations. However, BHF represents only the first- 
order approximation in the general hole-line expansion 
Conclusive studies 5, 6, 7] have shown that the hole- 
line expansion converges at third order (or second order 
with a continuous single-particle spectrum) and that such 
results are in good agreement with variational calcula- 
tions 3 of the binding energy per nucleon. Nonetheless, 
all such calculations have shown that it is very difhcult 
to obtain both the empirical saturation energy and sat- 
uration Fermi momentum simultaneously. In fact, such 
calculations using various models of the nucleon-nucleon 
interaction result in a series of saturation points which 
actually lie along a band, often referred to as the Coester 
band 3] , which deviates significantly from the empirical 
saturation point. For this reason it is now widely be- 
lieved that free-space two-nucleon interactions alone are 
insufficient to describe the properties of nuclear systems 
close to saturation density and that accurate results can 
only be achieved by introducing higher-order effects, e.g. 
three-nuclcon forces [l3| or relativistic effects [ll|. 
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In the present work, wc shall carry out calculations of 
the nuclear binding energy for symmetric nuclear matter 
using a framework based on a combination of the re- 
cently developed low-momentum NN interaction Viow-k 
[T^ . [isl . [3, [l^, [l^ UTl and the ring-diagram method for 
nuclear matter of Song et al. [18], which is a model-space 
approach where the particle-particle hole-hole (pphh) 
ring diagrams for the potential energy of nuclear matter 
are summed to all orders. In previous studies a model 
space of size A ^ 3 fm~^ was used to obtain improved 
results compared with those from the BHF method. Such 
an improvement can be attributed to the following desir- 
able features in the ring diagram approach. First, the 
ground-state energy shift AEq in the BHF approach is 
given by just the lowest-order reaction matrix (G-matrix) 
diagram (corresponding to diagram (b) of Fig. [1] with 
the dashed vertex representing G). It does not include 
diagrams corresponding to the particle-hole excitations 
of the Fermi sea. Such excitations represent the effect 
of long-range correlations. In contrast, the pphh ring 
diagrams, such as diagrams (c) and (d) in Fig. [Tl are 
included to all orders in the ring-diagram approach. Sec- 
ondly, the single-particle (s.p.) spectrum used in the ring- 
diagram approach is different from that in early BHF cal- 
culations, where one typically employed a self-consistent 
s.p. spectrum for momenta k < kp and a free-particle 
spectrum otherwise. Thus the s.p. spectrum had a large 
artificial discontinuity at kp. The s.p. spectrum used 
in the ring diagram approach is a continuous one. The 
importance of using a continuous s.p. spectrum in nu- 
clear matter theory has been discussed and emphasized 
in Ref . 0, 0] ■ Within the above ring diagram framework, 
previous calculations [l^ using G-matrix effective inter- 
actions and A ~ 3 fm~^ have yielded saturated nuclear 
matter that is slightly overbound {E/A ~ —18 MeV) 
and that saturates at too high a density {kp ~ 1.6 fm^^) 
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compared to empirical data. These results are consistent, 
within theoretical errors, with calculations based on the 
third-order hole-line expansion and variational methods 
(see Refs. 




(b) (c) (d) 

FIG. 1: Diagrams included in the pphh ring-diagram summa- 
tion for the ground state energy shift of symmetric nuclear 
matter. Included are (a) self-energy insertions on the single- 
particle propagator and (b) general pphh correlations. 

In the past, the above ring-diagram approach [H, [l^ 
employed the G-matrix interaction which is energy de- 
pendent, meaning that the whole calculation must be 
done in a "self-consistent" way. The calculation would 
be greatly simplified if this energy dependence, and thus 
the self-consistency procedure, were removed. Such an 
improvement has occurred in the past several years with 
the development of a low-momentum NN interaction, 
Viow-k, constructed from renormalization group tech- 
niques [H, [H [13, [H [H, [i3- As discussed in these 
references, the Viow-k interaction has a number of desir- 
able properties, such as being nearly unique as well as 
being a smooth potential suitable for perturbative many 
body calculations [2^ . Furthermore, Viow-k is energy in- 
dependent, making it a convenient choice for the interac- 
tion used in ring diagram calculations of nuclear matter. 

The Viow-k interaction has been extensively used in 
nuclear shell-model calculations for nuclei with a few va- 
lence nucleons outside a closed shell. As reviewed re- 
cently by Coraggio et al. the results obtained from 
such shell-model calculations are in very good agreement 
with experiments. However, applications of the Viow-k 
interaction to nuclear matter have been relatively few 
[21I [2^ m, [23|. A main purpose of the present work 
is to study the suitability of describing symmetric nu- 
clear matter using Viow-k- A concern about such ap- 
plications is that the use of Vioiv-k alone may not pro- 
vide satisfactory nuclear saturation. As illustrated in ref. 
[2l| . Hartree-Fock (HF) calculations of nuclear matter 
using Viow-k with a cutoff momentum of A ^ 2.0 fm~^ 
do not yield nuclear saturation-thc calculated E/A de- 
creases monotonically with kp up to the decimation scale 
A. 

In this work, we carry out a ring-diagram calculation of 
symmetric nuclear matter with Viow-k- We shall show in 
detail that satisfactory results for the saturation energy 
and saturation Fermi momentum can be obtained when 
one takes into account the following two factors: a suit- 



able choice of the cutoff momentum and the in-medium 
modification of meson masses. As we shall discuss, ring- 
diagram correlations at intermediate momenta (fc ~ 2.0 
fm~^) have strong medium dependence and are impor- 
tant for nuclear saturation. To include their effects one 
needs to use a sufficiently large decimation scale A so that 
the above momentum range is not integrated out. We 
have carried out ring-diagram calculations for symmetric 
nuclear matter using A ~ 3 fm^^ with several modern 
high-precision NN potentials, and the results yield nu- 
clear saturation. However, E/A and kp at saturation 
are both considerably larger in magnitude than the cor- 
responding empirical values. Great improvement can be 
obtained when one takes into account medium modifi- 
cations to the exchanged mesons. Clearly mesons in a 
nuclear medium and those in free space are different: the 
former are "dressed" while the latter are "bare" . Brown 
and Rho have suggested that the dependence of meson 
masses on nuclear density can be described by a simple 
equation known as Brown-Rho scaling [sil . [s^ : 
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where 17^1 is the axial coupling constant, $ is a func- 
tion of the nuclear density n, and the star indicates in- 
medium values of the given quantities. At saturation 
density $(no) ~ 0.8. In a high-density medium such 
as nuclear matter, these medium modifications of meson 
masses are significant and can render Vnn quite different 
from that in free space. Thus, in contrast to shell model 
calculations for nuclei with only a few valence particles, 
for nuclear matter calculations it may be necessary to use 
a Vjv N with medium modifications built in. In the present 
work, we shall carry out such a ring-diagram summation 
using a Brown-Rho scaled NN interaction. 

The Skyrme [34| interaction is one of the most success- 
ful effective nuclear potentials. An important component 
of this interaction is a zero-range 3-body force, which is 
equivalent to a density-dependent 2-body force. Note 
that the importance of 3-body interactions in achieving 
nuclear saturation with low-momentum interactions has 
been extensively discussed in the literature (see ref. (2^ 
and references quoted therein). In the last part of our 
work, we shall study whether the density dependence 
from Brown-Rho scaling can be well represented by that 
from an empirical density-dependent force of the Skyrme 
type. 

The organization of this paper is as follows. In Sections 
im and imi we outline our model space pphh ring-diagram 
calculation for the nuclear binding energy and the con- 
cept of Brown-Rho scaling respectively. In Sec. IIVI we 
present our computational results. A brief conclusion 
can be found in Sec. |Vl 



3 



II. SUMMATION OF pphh RING DIAGRAMS 

In this section we describe how to calculate the prop- 
erties of symmetric matter using the low-momentum ring 
diagram method. We employ a momentum model space 
where all nucleons have momenta A; < A. By integrating 
out the k > A components, the low-momentum inter- 
action Viow-k is constructed for summing the pphh ring 
diagrams within the model space. 

The ground state energy shift AEq = Eq — Eq'^'^ for 
nuclear matter is defined as the difference between the 
true ground-state energy Eq and the corresponding quan- 
tity for the non-interacting system Eq'^°. In the present 
work, we consider AEq as given by the all-order sum of 
the pphh ring diagrams as shown in (b), (c) and (d) of 
Fig. [11 

We shall calculate the all-order sum, denoted as AEq^, 
of such diagrams. Each vertex in a ring diagram is the 
renormalized effective interaction Viow-k corresponding 
to the model space k < A. It is obtained from the follow- 
ing T-matrix equivalence method [H, [H, [H, [H [H, ■ 
Let us start with the T-matrix equation 



T{k',k,k^) = V{k',k) 



(2) 



where is a bare NN potential. In the present work 
we shall use the CD-Bonn [13], Nijmegen-I [111 and 
Idaho(chiral) [29j NN potentials. Notice that in the 
above equation the intermediate state momentum q is 
integrated from to cxd. We then define an effective low- 
momentum T-matrix by 



Tion,-k{p' ,P,P^) = Viou,-k{p' ,p) 



^ Viow-k{p' , q)Tiow-kiq,p,p'^) 

— 



V / q^dq 
Jo 



(3) 



where the intermediate state momentum is integrated 
from to A, the momentum space cutoff. The low mo- 
mentum interaction Viow-k is then obtained from the 
above equations by requiring the T-matrix equivalence 
condition to hold, namely 

T(p',p,p2) = To»^fc(p',P,p'); (p',p) < A. (4) 

The iteration method of Lee-Suzuki-Andreozzi [13, [H, 
[26j has been used in obtaining the above Viow-k- 

With Viow-k, our ring diagram calculations are rela- 
tively simple, compared to the G-matrix calculations of 
ref. 18]. Within the model space, we use the Hartree- 
Fock s.p. spectrum calculated with the Viow-k interac- 
tion, and outside the model space we use the free particle 
spectrum. In other words, 



h^e/2m + Y.h<kA^h\Vion,~k\kh)- k<A 
h^k^/2m; k>A. 



(5) 



The above s.p. spectrum is medium (kp) dependent. 
Our next step is to solve the model space RPA equation 

^l(e.i + f-])kj,ef + Xifiifij ~ nin.j){ij\Viow-k\ef)] 

xy„(e/,A) =c^„r„(y,A); (z,i,e,/)<A, (6) 

where Ua — 1 for a < kp and na — for k > kp; also 
ha = {1 — na). The strength parameter A is introduced 
for calculational convenience and varies between and 1. 
Note that the above equation is within the model space 
as indicated by (i, j, e, /) < A. The transition amplitudes 
Y of the above equation can be classified into two types, 
one dominated by hole-hole and the other by particle- 
particle components. We use only the former, denoted 
by Ym, for the calculation of the all-order sum of the 
pphh ring diagrams. This sum is given by [isl. [23. Isol] 



m ijkKA 

xY;^ikl,X){ij\Vio,u-k\kl), (7) 

where the normalization condition for Ym is 
{Yjn\^\Y^) = -1 and Q{i,j) = {niUj - TiiUj). In 
the above, E means we sum over only those solutions of 

m 

the RPA equation ^ which are dominated by hole-hole 
components as indicated by the normalization condition. 

The all-order sum of the pphh ring diagrams as indi- 
cated by diagrams (b-d) of Fig. [T] is given by the above 
AEq^. Since we use the HF s.p. spectrum, each propa- 
gator of the diagrams contains the HF insertions to all 
orders as indicated by part (a) of the figure. Clearly our 
ring diagrams are medium dependent; their s.p. prop- 
agators have all-order HF insertions which are medium 
dependent, as is the occupation factor {hihj — niUj) of 
the RPA equation. 



III. BROWN-RHO SCALING AND IN-MEDIUM 
NN INTERACTIONS 

Nucleon-nucleon interactions are mediated by meson 
exchange, and clearly the in-medium modification of me- 
son masses is important for NN interactions. These mod- 
ifications could arise from the partial restoration of chi- 
ral symmetry at finite density/temperature or from tra- 
ditional many-body effects. Particularly important are 
the vector mesons, for which there is now evidence from 
both theory [H, Ha, [13] and experiment [H, js^ that the 
masses may decrease by approximately 10 — 15% at nor- 
mal nuclear matter density and zero temperature. This 
in-medium decrease of meson masses is often referred to 
as Brown- Rho scaling [3l|, [12] . For densities below that 
of nuclear matter, it is suggested [1^ that the masses 
decrease linearly with the density n: 



my 
mv 
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(8) 



4 



where rriy is the vector meson mass in-medium, ng is 
nuclear matter saturation density and C is a constant of 
value -0.10-0.15. 

We study the consequences for nuclear many-body 
calculations by replacing the NN interaction in free 
space with a density-dependent interaction with medium- 
modified meson exchange. A simple way to obtain such 
potentials is by modifying the meson masses and relevant 
parameters of the one-boson-exchange NN potentials 
(e.g. the Bonn and Nijmegen interactions). The satura- 
tion of nuclear matter is an appropriate phenomenon for 
studying the effects of dropping masses [H, [s^ , since the 
density of nuclear matter is constant and large enough 
to significantly affect the nuclear interaction through the 
modified meson masses. 

One unambiguous prediction of Brown-Rho scaling in 
dense nuclear matter is the decreasing of the tensor force 
component of the nuclear interaction. The two most im- 
portant contributions to the tensor force come from tt 
and p-meson exchange, which act opposite to each other: 

Vffir) = -^"^pTl • T2Si2f3impr), (9) 
Vj{r) = ^m^r, ■ T2Si2h{m^r), (10) 



with Tav = \{ri+r2). 

The general structure of Vskyrme is rather similar to 
the effective interactions based on effective field theo- 
ries (EFT) [ill, '^itt' ^{hi) corresponding to Viow-k and 
V{i,j,k) to the EFT 3-body force. The Skyrme 3-body 
force, however, is much simpler than that in EFT. We 
shall compare in the next section the density dependent 
effect generated by the medium modified NN interaction 
with that from an empirical 3-body force of the Skyrme 
type. 



IV. RESULTS AND DISCUSSIONS 

In this section, we shall report computational results 
for the binding energy of symmetric nuclear matter cal- 
culated with an all-order summation of low-momentum 
pphh ring diagrams. The method is already outlined and 
discussed in the above sections. As mentioned above, 
we employ a model space approach. Starting from var- 
ious bare NN interactions, we first construct the low- 
momentum interactions Viow-k with a particular choice 
of the cutoff momentum A. The low-momentum (< A) 
pphh ring-diagrams are then summed to all-orders as 
given by Eq. ([7]) to give the binding energy. 



/3(mr) = + ^ + e— . (11) 

In Brown-Rho scaling the p meson is expected to decrease 
in mass at finite density while the pion mass remains 
nearly unchanged due to chiral invariance. Therefore, 
the overall strength of the tensor force at finite density 
will be significantly smaller than that in free space. As 
we shall discuss later, this decrease in the tensor force 
plays an important role for nuclear saturation. 

The Skyrme effective interaction has been widely used 
in nuclear physics and has been very successful in de- 
scribing the properties of finite nuclei as well as nuclear 
matter |34l. This interaction has both 2-body and 3-body 
terms, having the form 

i<j i<j<k 

Here V{i,j) is a momentum (fc) dependent zero-range in- 
teraction, containing two types of terms: one with no mo- 
mentum dependence and the other depending quadrati- 
cally on k. V{i,j) corresponds to a low-momentum ex- 
pansion of an underlying NN interaction. Its 3-body term 
is a zero-range interaction 

Vit.j, k) = t:i5{n - fj)5(fj - Tk) (13) 

which is equivalent to a density-dependent 2-body inter- 
action of the form 

Vp{l, 2) = \h5{Ti - r2)p{rav) (14) 



A. Single-particle spectrum and nuclear binding 
energy 

First, we shall look carefully into the role of A in our 
ring-diagram calculation. Let us start with the single 
particle (s.p.) energy e/c. Obtaining is the first step 
in our ring-diagram calculation. Within our model space 
approach, tk is given by the Hartree-Fock spectrum for 
fc < A, while for k > A, ek is taken as the free spectrum 
(see Eq. ([S])). As emphasized before, the s.p. spectrum 
obtained in this way will in general have a discontinuity 
at A. Such a discontinuity is a direct consequence of 
having a finite model space. It is of much interest to 
study the s.p. spectrum as A is varied. In Fig. ^ we 
plot the spectrum for different values of A ranging from 
2-4 fm~^ We observed that with A = 2.0 fm"^, the 
discontinuity at A is relatively large; there is a gap of 
about 50 MeV between the s.p. spectrum just inside A 
and that outside. However, this discontinuity decreases 
if A is increased to around 3 fm~^. At this point, the 
s.p. spectrum is most "satisfactory" in the sense of being 
almost continuous. A further increase in A will result in 
an "unreasonable" situation where the s.p. spectrum just 
inside A becomes significantly higher than that outside. 
This is clearly shown in the data of A = 4.0 fm~^. The 
above results suggest that to have a nearly continuous s.p. 
spectrum, which is physically desirable, it is necessary to 
use A — 3 fm~^. 

Next, we shall look into the effect of A on the nuclear 
binding energy. Once the s.p. energies are obtained, the 
all-order ring-diagram summation can be carried out (see 
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FIG. 2: Dependence of the model-space s.p. spectrum on the 
decimation scale A for symmetric nuclear matter at the em- 
pirical saturation density. The CD-Bonn potential is used in 
the construction of Viow-k- 

Eqs. ^ and ([7])). Let us first discuss the computational 
results based on the CD-Bonn potential. Results from 
various A ranging from 2 — 3.2 fm^^ are shown in Fig. 
121 Let us focus on (i) the overall saturation phenomena 
and (ii) the numerical values of the binding energy and 
the saturation momentum. 

(i) We observe that the nuclear binding energy exhibits 
saturation only when A is ~ 3 fm~^ and beyond. This 
reflects the importance of ring diagrams in the intermedi- 
ate momentum region (fc ~ 2 fm^^). To illustrate, let us 
compare the results for the cases of A = 2 and 3 fm"-'^ . As 
indicated by Eqs. (QHll), Viow-k includes only the k > A 
pp ladder interactions between a pair of "free" nucleons; 
there is no medium correction included. Thus the above 
two cases treat correlations in the momentum region be- 
tween 2 and 3 fm"-'^ differently: the former includes for 
this momentum region only pp ladder interactions with 
medium effect neglected, while the latter includes both 
pp and hh correlations with medium effect, such as that 
from the Pauli blocking, included. Our results indicate 
that the medium effect in the above momentum region is 
vital for saturation. 

For nuclear matter binding energy calculations, there 
is no first-order contribution from the tensor force (Vr); 
its leading contribution is second order of the form 
{^Si\Vt^Vt\^Si) where Q stands for the Pauh blocking 
operator and e the energy denominator. Thus the con- 
tribution from the tensor force depends largely on the 
availability of the intermediate states; this contribution 
is large for low kp but is suppressed for high kp- To illus- 
trate this point, we plot the potential energy of nuclear 
matter from the ^Sq and ^Si — '^Di channels separately 
in Fig. [4] The behavior of the potential energy in these 
two channels differ in a significant way. The ^Sq channel 
is practically independent of the choice of A, as displayed 
in the upper panel of the figure. This indicates that for 



this channel the effects from medium corrections and hh 
correlations are not important. Also the PE/A from this 
channel does not exhibit saturation at a reasonable kp- 
In the lower panel of the figure, we display the PE/A for 
the '^iS'i— '^Di channel where the tensor force is important. 
As seen, PE/A does not exhibit saturation when using 
A = 2 fm~^. On the contrary, the result using A = 3 
fm^^ shows a clear saturation behavior. This is mainly 
because that in the former case the Pauli blocking effect 
is ignored for the momentum region 2 — 3 fm~^ while it 
is included for the latter. To have saturation, we should 
not integrate out the momentum components in the NN 
interaction that are crucial for saturation. Considering 
also the effect of A on the s.p. spectrum, we believe that 
A = 3.0 fm~^ is a suitable choice for our ring-diagram 
nuclear matter calculation. Notice that a model space 
~ 3 fm^^ has been used in other similar ring summation 
calculations using G- matrix effective interactionfisl [l9j . 
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FIG. 3: Results for the energy per nucleon (E/A) of sym- 
metric nuclear matter obtained by summing up the pphh ring 
diagrams to all orders. Low-momentum NN interactions, con- 
structed from the CD-Bonn potential, with various cutoffs A 
are used in the ring diagram summation. 

(ii) We have performed a similar ring summation with 
the Nijmegen I and Idaho potentials. Results with 
A — 3.0 fm"^ are compared with that from CD-Bonn as 
shown in Fig. [5l The saturation energies for these three 
potentials are located between —19 and —23 MeV, while 
the saturation momentum ranges from 1.75 to 1.85 fm~^. 
These quantities are considerably larger than the empir- 
ical values of —16 MeV and 1.4 fm^^, respectively. We 
believe that improvements can be obtained if one takes 
into account the medium dependence of the NN inter- 
action. Namely, instead of using a Viow~k constructed 
from a bare NN interaction, one should employ a Viow-k 
constructed from a "scaled" NN interaction according to 
the nuclear density. Below we shall report how we incor- 
porate such effects into our ring diagram summation. 
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FIG. 4: Potential energy per nucleon (PE/A) in the ^So and 
^Si — ^Di channels of symmetric nuclear matter from sum- 
ming up pphh ring diagrams to all orders. The CD-Bonn 
potential is used in the construction of Viom-k- 



Brueckner-Hartree-Fock formalism showed that there is 
too much attraction when the a meson is scaled according 
to dH). However, a microscopic treatment of a me- 
son exchange in terms of correlated 27r exchange showed 
that the medium effects on the a are much weaker than 
in ([8]). Therefore, in our ring diagram summation using 
the Brown-Rho scaled Nijmegen II interaction, we em- 
ploy a range of scaling parameters Co- between 0.075 and 
0.09. Our calculations are shown in Fig. [6l With Brown- 
Rho scaling, the numerical values for both the saturation 
energy and saturation momentum are greatly improved. 
Whereas the unsealed potential gives a binding energy 
BE /A ~ 20 MeV and A:^ ~ 1.8 fm^^, the scaled poten- 
tial gives BE/A ~ 14 - 17 McV and k% ~ 1.30 - 1.45 
fm~^ for a a meson scaling constant Co- ^ 0.08 — 0.09, 
in very good agreement with the empirical values. We 
conclude first, that the medium dependence of nuclear 
interactions is crucial for a satisfactory description of nu- 
clear saturation and second, that within the framework of 
one-boson-exchange NN interaction models one can ob- 
tain an adequate description of nuclear matter saturation 
by including Brown-Rho scaled meson masses. 




1 1.25 1.5 1.75^ 2 2.25 2.5 8 1 1.2 1.4 1.6 1.8 2 

[fm'l [fm'] 



FIG. 5: The binding energy of symmetric nuclear matter 
from the low-momentum ring-diagram summation using var- 
ious NN potentials. A momentum-space cutoff of A = 3.0 
fm~^ is used. 



B. Nuclear binding energy with Brown-Rho scaling 

The concept of Brown-Rho scaling has already been 
discussed in Sec. IIIIl The medium effects on the NN in- 
teraction resulting from the in-medium modification of 
meson masses shall have a profound effect on nuclear 
binding. To incorporate this in our ring-diagram cal- 
culation we work with the Nijmegen potential, which is 
one of the pure one-boson-exchange NN potentials. The 
bare Nijmegen is first Brown-Rho scaled (see Eq. ([8])) 
with the dropping mass ratio C chosen to be 0.15. Vec- 
tor meson masses in a nuclear medium have been widely 
studied both theoretically and experimentally, but the a 
meson mass is not well constrained. Previous calcula- 
tions I33l of nuclear matter saturation within the Dirac- 



FIG. 6: (color online). The binding energy of symmetric nu- 
clear matter from the Brown-Rho scaled low-momentum Ni- 
jmegen II interaction using the ring-diagram summation with 
A = 3.0/m~^. Calculations for different choices of the a me- 
son scaling constant C'a are shown. 



C. Nuclear binding energy with 3-body force of 
the Skyrme type 

As discussed earlier in section III, the widely used 
Skyrme interaction contains a 3-body term that is equiv- 
alent to a density-dependent 2-body interaction. It is of 
much interest to study whether our result with Brown- 
Rho scaled Nijmegen potential can be reproduced with 
the unsealed Nijmegen plus an effective 3-body inter- 
action of the Skyrme type which is characterized by a 
strength parameter ^3 (see Eq. (14)). In Fig. [7] we com- 
pare the results using ^3 = 1250 with our previous calcu- 
lations using the Brown-Rho scaled Nijmegen II potential 
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with a a meson scaling constant of Ca = 0.087. In all 
calculations A = 3.0 is used. We note that satisfactory 
results for the saturation energy and Fermi momentum 
are obtained using either Brown-Rho scaling or a 3NF of 
the Skyrme type. However, the nuclear incompressibility 
is considerably larger in the case of Brown-Rho scaling. 
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FIG. 7: The binding energy of symmetric nuclear matter from 
the low momentum ring diagram summation with A = 3.0 
fm~^. Three different interactions are used: (1) the medium- 
independent Nijmegen II interaction, the Brown-Rho scaled 
interaction with = 0.087, and finally the Nijmegen II in- 
teraction supplemented with a 3NF of the Skyrme type with 
U = 1250. 



hole (pphh) ring diagrams within a momentum model 
space fc < A were summed to all orders. The signifi- 
cant role of the intermediate momentum range ('^ 2.0 
fm~^) for nuclear saturation was discussed. We con- 
cluded that in the ring diagram summation, having a 
sufhciently large model space is important to capture 
the saturation effect from the intermediate momentum 
components. Various bare NN potentials including CD- 
Bonn, Nijmegen and Idaho have been employed, resulting 
in nuclear saturation with A = 3.0 fm~^. However, the 
resulting binding energy and saturation momentum are 
still much larger than empirical values. Improvement can 
be obtained when we take into account the medium mod- 
ification of NN interaction. We first constructed Viow-k 
from a medium-dependent Brown-Rho scaled NN poten- 
tial and then implemented this into the ring-diagram 
summation. Satisfactory results of E/A ~ —15 MeV 

and kp'' ~ 1.4 fm~^ could then be obtained. We showed 
that these saturation properties are well reproduced by 
the first ring-diagram approach with the addition of an 
empirical 3-body force of the Skyrme type. 

In the future, it is of much interest to carry out a BCS 
calculation on nuclear matter with Viow-k, particularly 
for the '^S'l —^Di channel where earlier calculations using 
bare NN interactions revealed a gap of 10 MeV around 
normal nuclear matter densities [40, [4l|. Recently, Viow-k 
has been applied to obtain the equation of state of neu- 
tron matter [Ulill and the '^Sq pairing gapSlHl]. A 
similar calculation on nuclear matter which incorporates 
the tensor correlations is obviously important and we 
plan to investigate it in the future. 



V. CONCLUSION 

We have studied the equation of state for symmet- 
ric nuclear matter using the low-momentum nucleon- 
nucleon (NN) interaction Viow-k- Particle-particle hole- 
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